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We calculate the ground state and simulate the dynamics of a finite chain of spins with Ising 
nearest-neighbor interactions and a Dicke collective spin interaction with a single mode cavity field. 
We recover the signatures of first and second order phase transitions predicted by mean field theory, 
while for small chains, we find significant and non-trivial finite size effects. Below the first order 
phase transition, even quite large spin chains of 30-40 spins give rise to a mean photon number 
and number fluctuations significantly above the mean field vacuum result. Near the second order 
phase critical point, our calculations reveal photon number fluctuations that grow beyond Poisson 
statistics with the size of the spin chain. We simulate the stochastic evolution of the system when 
the cavity output field is subject to homodyne detection. For an initial state close to the first order 
phase-transition the random character of the measurement process causes a measurement-induced 
symmetry-breaking in the system. This symmetry-breaking occurs on the time-scale needed for 
an observer to gather sufficient information to distinguish between the two possible (mean-field) 
symmetry-broken states. 
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I. INTRODUCTION 

The problem of a single quantum system coupled to 
a large number of mutually interacting systems is found 
in a wide range of physics problems: an electron spin 
in a quantum dot or in an NV-center interacting with 
nuclear spins, an atom interacting with a continuum of 
quantized field modes, a cavity field mode interacting 
with an ensemble of absorbing atoms. 

In this manuscript we present a matrix product state 
(MPS) analysis of the quantum states of a many body 
system, combining a one dimensional chain of spin-1/2 
particles with nearest neighbor interaction and a cavity 
oscillator mode with a quantized radiation field which 
is linearly coupled to all the spins. These two interac- 
tions are separately known as the Ising and the Dicke 
models with known many-body properties and phase di- 
agrams, and models which combine the spin-spin and 
spin-oscillator models interactions have been studied by 
mean field theory in recent publications [1-3]. By use 
of the matrix product state analysis we will be able to 
study non-mean field and finite size effects in such highly 
relevant theoretical models. 

Numerous examples of quantum systems with local in- 
teraction in one spatial dimension have revealed a de- 
crease of entanglement with the distance between pairs 
of particles. Inspired to a large extent by insights from 
quantum information science, matrix product states have 
been developed as a very powerful method for the rep- 
resentation and solution of quantum many-body prob- 
lems in physics. The basic idea is to explicitly use the 
restriction of the entanglement between remote compo- 
nents of the full quantum system in a suitable variational 
ansatz for the quantum state. A review of the devel- 
opment of these ideas, examples of their application to 
different systems, and their equivalence to the so-called 
density matrix renormalization group (DMRG) is offered 



in [4, 5]. The interaction of the the cavity field with all 
the other subsystem degrees of freedom seemingly vio- 
lates the usual requirement for matrix product state to 
apply. We show, however, that by an adequate represen- 
tation we can indeed implement and validate the method 
even with a globally coupled mode. 

The paper is organized as follows. In Sec. II, we re- 
call the basic idea of the matrix product state method, 
and we describe how one may include a non-local quan- 
tum degree of freedom, and why this does not violate 
essential properties of the method. In Sec. Ill, we in- 
troduce the Dicke-Ising Hamiltonian, and we present our 
numerical results for the ground state of the system as 
a function of the interaction parameters. In symmetry 
breaking mean field theory, first and second order phase 
transitions have been identified in the model [1, 2], and 
we show how these transitions reveal themselves in sys- 
tems of finite size without breaking of the symmetries of 
the Hamiltonian. Finally, in Sec. IV, we show how mea- 
surements on the oscillator mode can be simulated within 
the matrix product state method, and how the random 
measurement outcomes in a natural manner drive the 
system into symmetry broken states. In Sec. V, we sum- 
marize and conclude. 



II. THE MATRIX PRODUCT STATE METHOD 

For one-dimensional lattice systems, matrix product 
states (MPS) [4-7] provide an efficient parametriza- 
tion of the many-body Hilbert space. For a system 
with L sites and a conventional tensor product basis 
® 1*2) ® ■ ■ ■ \ih) = |*i ■ • • *l) the matrix product state 
parametrization of the expansion coefficients takes the 
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For each basis state \ij) on a site j one defines a, D x D 
matrix A^' lj , or, equivalently, for each site j one de- 
fines a D x D x d tensor, , where d is the single-site 
Hilbert space dimension. By multiplying these matrices 
together and extracting the trace as specified in (1), we 
thus obtain a set of expansion coefficients on the ten- 
sor product basis states \ii . . .i^). The total number of 
parameters in all matrices equals LD 2 d and grows only 
linearly with the number of sites, unlike the number of 
tensor product basis states which grows as d L . The valid- 
ity of the expansion thus relies strongly on the restriction 
of the quantum system to a small part of its Hilbert space 
caused, in parts, by the locality of interactions. 

For D = 1 all expansions coefficients factor into prod- 
ucts of single site amplitudes, A^' lj , and the state in 
(1) is simply a product state, valid for non-interacting 
and non-entangled particles. For higher values of the 
matrix dimension D, however, the parametrization (f) 
forms the basis for a variational calculation of the sys- 
tem properties, which supports an increasing amount of 
entanglement between the subsystems. There is a host of 
theoretical work regarding which states matrix product 
states can represent efficiently (i.e. with a reasonably 
low D) [8-12]. In general the results reveal that one- 
dimensional systems with nearest-neighbor interactions 
and, in particular, gapped one-dimensional systems can 
be represented efficiently. It is also possible, however, to 
efficiently describe systems with interactions with both 
exponential and polynomial dependence on the distance 
between the particles [13, 14]. A key advantage of ma- 
trix product states is the ability to evaluate the energy 
variance (ip\H 2 \ip) - (if>\H\ij) 2 which gives an absolute 
bound on the accuracy of the state-energy, and also a very 
reliable indicator for convergence even for non-standard 
problems. 

In this paper we investigate a more general system than 
a one-dimensional chain. We will introduce one more de- 
gree of freedom which is allowed to interact simultane- 
ously with all the sites of the chain. This situation occurs 
in e.g. the Dicke-model [2, 15, 16] and Spin-star models 
[17-20]. The added degree of freedom mediates a cou- 
pling and hence permits correlations and entanglement 
between remote spins which is independent of their sep- 
aration. This, however, should not discourage use of the 
MPS method, which indeed does not rely on the particu- 
lar fall-off of correlations with distance but rather on the 
entanglement being finite between any two parts, e.g., 
defined by a cut at any given location in the chain. The 
fact that the global mode mediates interactions between 
any such parts of the spin chain may not cause a prob- 
lem, as long as their entanglement is kept finite. Our 
study, in parts, aims to investigate how crucial are the 
resulting limitations on the Hilbert space dimension and 



the degree of excitation of the global system. 

The natural way to represent our system is to de- 
compose the system Hilbert space into a spin chain-part 
and the external, global part W = %i a ttice ® %ext ■ The 
spin-chain space is a tensor product of L copies of a d- 
dimensional space T-L lattice = %2te- The system state can 
then be expanded as a spin chain-state for each basis- 
state of the global degree of freedom, e.g. 
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where |n) is the n'th basis-state of the external degree 
of freedom and \ip n ) could then be an MPS with tensors 
Al*,n] defined for each site i and each global system basis 
state n. If the external degree of freedom is a a spin 
of magnitude S, d cxt = 25 + 1, while for a harmonic 
oscillator a truncation would be necessary and d C xt will 
then denote the cut-off excitation number. 

By choosing \tp n ) as MPS-states we can also rewrite 
the full state \ip) as a block-diagonal MPS-state, i.e. 

|V) = J2 ■&(" n ^ I1],<1 ■■■A [L] ' iL ) \n,i u ...i L ) (3) 
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and v n is a vector of zero matrices except the block cor- 
responding to the n'th component of A^ which is an 
identity. Note that (3) is equivalent to a standard MPS 
ansatz, only it is restricted to a certain block-diagonal 
structure, with a new effective D = Dd ext . In practice, 
however, we do not need to restrict to this block-diagonal 
form. By allowing general matrices we can have a more 
compact representation when, e.g. the two subsystems 
are in a product state. 

For systems with a few degrees of freedom interacting 
with a spin chain the Hamiltonian will most naturally 
contain terms which are local in the chain and the ex- 
ternal degrees of freedom, or which are simple products 
of such operators. All such terms have compact matrix 
product operator expressions and the tensor networks 
arising in matrix-elements can be evaluated efficiently. 
For example, an interaction term V = A J2i=i where 
A acts on "H ext and Xi acts on the i'th lattice site can 
be written as a matrix product operator (MPO) with 
tensors 

M° = (A) M 1 = (X, I) , W = ( x. J) . mL = (*/) 

where A is a d ex t x d ex t matrix, and the Xi are d x d ma- 
trices. These compact representations of operators are 
key to the application of the variational methods devel- 
oped for matrix product states; both for eigenstates [4] 
and time-evolution [21]. 
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For unitary time-evolution we will use the Krylov 
propagation algorithm [21] instead of the more common 
Suzuki- Trotter decomposition. The Krylov method con- 
sists of constructing a basis for the subspace spanned 
by , H \ip) , H 2 |</>) , • • • , H N |V>)}- By truncating at a 
finite power of the Hamiltonian we can get a good ap- 
proximation to the matrix exponential of H and thus 
perform time-evolution very efficiently within this sub- 
space. Once the state has been propagated within 
the Krylov subspace, it is injected into the full Hilbcrt 
space again. 

In order to perform these operations for an MPS we 
need to be able to construct the Krylov basis and inject 
the propagated vector into the Hilbert space using varia- 
tional methods. There is an efficient algorithm for finding 
an MPS of a given virtual dimension which best approx- 
imates an operator-MPS product, i.e. find |y) such that 
j| |x) — T \ip)\\ 2 is smallest. This works essentially as long 
as the expression (x\T\ip) can be evaluated efficiently. 
The orthogonal Krylov basis vectors arc given by the 
Gram-Schmidt algorithm as 

i=0 x 1 ' 

for j = 1, . . . N — 1 and |0) = \ip). Since we can contract 
terms of the form (i\H\j) when and \j) are matrix 
product states efficiently, we can also find an MPS which 
is closest to each of the Krylov basis vectors efficiently. 
After performing the unitary time evolution within the 
Krylov subspace we have a state U\ip) w Y^^o c j \ j) 
for which the nearest MPS is found using the the same 
variational technique. 

The ground state calculations presented in this 
manuscript were all performed with sweeps of local min- 
imization of the effective Hamiltonian [4] (using D <~ 
128), whereas the time-evolution was calculated using 
the Krylov subspace technique described in [21] (using 
D~512). 



III. THE DICKE-ISING MODEL 

We describe our interacting spin chain and the global 
oscillator by the Hamiltonian, 

L L-1 L 

(4) 

where cr| , £ = x,y,z denote Pauli spin-1/2 operators for 
the j th spin, w,a,a^ denote the eigenfrequency, and the 
lowering and raising operators for the oscillator, J is the 
Ising coupling constant, g/VL is the Dicke coupling con- 
stant, normalized to the system size, and ft. is a transverse 
bias (magnetic) field. The goal of our calculations is to 



identify the phase diagram under variation of the control 
parameters J,h,g in units where u> = 1. 

The analyses in [1, 2] suggest that in the thermody- 
namic limit L — > oo a mean field description of the os- 
cillator mode should be valid. The mean field approxi- 
mation essentially implies that the oscillator-spin chain 
system state is factored into a coherent state \a) and an 
Ising-interacting chain with an effective transverse mag- 
netic field h 2 s = h 2 + 4g 2 $i(a) 2 . We expect that the 
MPS ansatz (3) should work well in this limit, since it 
only requires tensor-dimension 1 on the link connecting 
the oscillator and the spin chain. 

The Hamiltonian (4) has a discrete symmetry P = 
it ® o\ . . . (7% where n^cnr = —a. This symmetry is bro- 
ken in the mean field description since 7r \ a) = \—a) and 
the free energy in the thermodynamic limit is invariant 
under the replacement a — >• —a. For a finite-size sys- 
tem the symmetry is not broken and the system ground 
state is either a +1 or —1 cigenstate of P. Since the 
mean field theory predicts a proportional to V~L, and 
(a\ — a) = cxp(— 2 |a| 2 ), where matrix elements between 
| a) and |— a) of a^a and a + a^ will be exponentially sup- 
pressed with increasing L, it becomes exceedingly dif- 
ficult, even for moderate systems sizes, to numerically 
identify states with the proper symmetry of the Hamil- 
tonian without explicit attention to symmetry. 

We have used the MPS method to calculate the prop- 
erties of the lowest eigenstates of the Hamiltonian (4) 
for L = 32 spins. In Figure 1 we have plotted the val- 
ues of different physical quantities in four simplex plots. 
The simplex plots restricts the four control parameters 
by fixing one of them, e.g., us = 1, and by imposing a 
linear constraint among the others, such as h + g + J = 1 
(assuming also that h,g,J > 0). The parameter depen- 
dence of the given physical quantities are then shown by 
color (hue) in the resulting 2D plot. In any corner of 
the plots, one parameter attains its maximum and the 
others vanish, while along the opposing edge the same 
parameter vanishes while the other two vary under the 
constraint that their sum is fixed. Along the ft-J-cdge 
(the vertical edge in Figure 1) the coupling between the 
oscillator and the spins vanishes, and the model reduces 
to the conventional Ising model, while for the lower (h-g) 
edge, the model has no spin-spin interactions and there- 
fore corresponds to the Dicke phase transition model. 

The first panel in Figure 1 shows the expected excita- 
tion of the oscillator divided by the number of spins, and 
we observe a strong similarity with the results of the low 
temperature mean field theory (see Figure 1 in [2]). Due 
to the finite size of the system, phase transitions will be 
smoothed, but still, we do recognize the characteristics of 
the Dicke second order phase transition from the normal 
to the super-radiant state for low spin-spin interaction 
strength. The similarly smoothed characteristics of the 
Ising-transition can also be seen in the fourth panel. The 
spins change from pointing parallel to the z-axis to hav- 
ing no net magnetization but strong y — y-corrclations, 
((T*(T^) ~ e - !' - - 7 '!/^. Finally, the first order transition pre- 
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Figure 1: (Color online) From left to right: expected photon density (n) /L, Energy gap to the excited state, Energy 
gap to the second excited state and the correlation- length for a y . Values are for 32 spins with oj = 1, and 

h+J+g= 1. 
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Figure 2: (Color online) Photon number mean (upper 
panel) and variance (lower panel) as functions of h for 

ui = 1, g = 0.4, J = 0.1. Note the peak in Varn/rt 
indicating a super-poissonian distribution at the phase 
transition. As a function of L the peak height scales as 

Var n/n oc L 019 and the peak location scales as 
\h — h c \ oc L~ 5 where h c is the location of the critical 
point. 



Figure 3: (Color online) w = l,g = 0.49, J = 0.3. The 

upper panel shows the mean photon number as a 
function of h near the first order transition. The lower 
panel shows the photon number variance in the region 

from h = 0.1 to h = 0.25. The dashed curve in the 
upper figure shows the mean photon number for a low, 
but finite, temperature. 



dieted by [1, 2] can also be seen along the upper edge in 
the simplex. 

The excited state energy gaps (panels two and three 
in Figure 1) in our finite system calculations also reflect 
the significant difference between the first and second or- 
der transitions. The second order transition occurs from 
the J < h region of the phase diagram and the system 
therefore enters from a unique ground state with a fi- 
nite total magnetization in the z-direction. During the 



second order transition the two lowest eigenstates trans- 
form from spins pointing in opposite z-directions to the 
two degenerate mean-field states. 

In Figure 2 we investigate the effects of the finite size 
of the system in more detail where the mean excitation 
of the oscillator, normalized to the number of spins and 
its variance are plotted as functions of external magnetic 
field h with uj = 1, g = 0.4 and J = 0.1 for a vari- 
ety of system sizes. According to the mean field theory, 
a second order phase transition into the super-radiant 
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state occurs at h rts 0.312 with the q = (a + a) ) /v^-field 
quadrature as order parameter. This is indicated with a 
full black line in the upper panel. 

For small system sizes the system gradually changes 
from supporting a few to a larger number of oscillator 
excitation quanta per particle and the MPS results grad- 
ually approach the ones of the mean field theory. The de- 
gree of excitation quickly approaches the mean field value 
as the system size increases and already for L = 32 spins 
the system behaves similarly to the mean field theory. 
An important difference between the mean field theory 
and the numerical calculations is the normalized variance 
Var n/n shown in the lower panel. The mean field theory 
predicts a Poisson distribution for the photon number 
for all parameter values implying Var n/n = 1; in the 
numerical calculations which go beyond the mean field 
description, there is a substantial deviation from this re- 
sult near the second order transition where we clearly see 
a peak in the normalized excitation variance scaling as 
Var n/n cx L 019 for this set of parameters. 

This peak is most likely due to the oscillator ground 
state gradually splitting into two coherent components 
giving rise to a widening of the excitation number dis- 
tribution. In the normal phase the cavity is close to a 
vacuum state and develops first into a weakly excited 
state consisting primarily of the symmetric combination 
| a) + |— a), with an excitation number variance which 
is higher than the individual coherent-state components. 
As the system enters the super-radiant state the field 
components separate and become correlated with orthog- 
onal states of the spin chain. The field state then be- 
comes equivalent to an even statistical mixture of the 
two symmetry-broken states with conventional Poisson 
statistics. 

In Figure 3 we have plotted the photon number den- 
sity and Var n/n near a first-order transition. Near the 
phase transition one can discern oscillations in the pho- 
ton number due to the ground state and first excited 
state alternating between two close-lying levels with dif- 
ferent symmetry. At finite temperature these oscillations 
are washed out as can be seen by following the dashed 
curves in Figure 3. 

Accompanying the oscillations in the mean photon 
number, significant deviations from the Poissonian statis- 
tics expected from the mean-field theory, can be seen in 
the lower part of Figure 3. A striking feature in this pa- 
rameter regime is the significant effect of the finite size of 
the system. Here, L > 40 is needed in order to see a sharp 
drop in photon number below the critical point, and the 
photon number variance is very far from the mean-field 
result for even larger systems. In addition, the location 
of the first-order transition predicted by mean-field the- 
ory is only approached very slowly as a function of the 
system size. 
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Figure 4: (Color online) Time evolution of system being 
monitored by homodyne detection, (a) The conditioned 
expectation value of the a + a) quadrature of the cavity 
field. The dashed line are the values of q for the two 
symmetry-broken states in the thermodynamic limit, 
(b) The thick line is the von Neuman entropy of the 
average time-evolved oscillator state, and the thin lines 
are the entropy of the oscillator state in a few sample 
trajectories, (c) The expectation-value of the symmetry 
operator P. 



IV. TIME EVOLUTION DURING HOMODYNE 
DETECTION 

While the Hamiltonian commutes with the symmetry 
operator P and the energy eigenstates can be chosen as 
eigenstates of P, the symmetry breaking ansatz associ- 
ated with either one of the coherent states is a signifi- 
cant simplification of the problem and it makes accurate 
predictions in the macroscopic limit. The spontaneous 
symmetry breaking ansatz may be justified by the small- 
ness of any perturbation needed to make the symmetry 
broken state the actual ground state of the system, but it 
may also be justified by the coupling and entanglement of 
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the system with its environment. A measurement taking 
place within the environment gives a random outcome 
by the laws of quantum mechanics, and this randomness 
may, indeed, be sufficient to break the symmetry of the 
quantum state. 

In our case, it is natural to imagine that the cavity is 
lossy, and thus the field may leak out of the cavity and be- 
come subject to a phase sensitive measurement probing. 
We consider a situation where the emitted cavity field is 
measured with a homodyne detection scheme. The sys- 
tem then evolves according to the stochastic Schrodingcr 
equation [22-24] 



-iHdt 



^ ^a^a — (a + a 1 ) a + ^ (a + a^) 2 ^ dt 
+VR(a- X - {a + o)))dW |V) 



where (a+a < )/\/2 is the cavity quadrature component be- 
ing measured by the homodyne detector, dW is a Wiener 
increment and k is the cavity loss rate. 

To simulate this equation, we treat the unitary and 
stochastic dynamics using a split-step approach as in 
[25]. The system is first propagated within a Krylov 
subspace generated by H for a time At and the stochas- 
tic part can be propagated by selecting a normally dis- 
tributed random number Ay with variance At, and mean 
y/K(a + a^)dt and update the MPS-state according to 
the action of the operator 

Sl(Ay) = 1 - -aladt + ^fnaAy 

followed by renormalization [22]. This mimics the mea- 
surement back action and samples the distribution of ob- 
served signals Ay{t). The measurement back action acts 
locally on the system, and hence the determination of the 
resulting MPS state is efficient [25]. 

In Figure 4 we show a small ensemble of simulated 
trajectories for a system of 32 spins, co = 1, h = 0.2, 
J = 0.3 and g — 0.49 and n — 0.5. In Figure 4(a) the thin 
lines show the conditioned time-evolution of the expected 
value of the cavity quadrature q = (a + a))/ \[2 along 
with the value estimated from the two symmetry-broken 
thermodynamic states (dashed line). The ^-quadrature 
starts out at zero, but quickly stabilizes to a value of 
about ±2.1. This is consistent with the interpretation 
that by probing the cavity field the observer is able to 
distinguish the two symmetry-broken states and thereby 
drives the system to the state consistent with the obser- 
vation. 

The time-scale on which the states are distinguished 
corresponds nicely with the time-scale for the decay of 
a mesoscopic superposition state [26]. This time-scale 
is given by r w l/K|Aa| 2 , where Aa is the seperation 
between the two coherent components. For k = 0.50a; 
and |Aa| w 4.2 we estimate ujt rts 0.11, which is con- 
sistent with the time-scale observed in Figure 4. The 
thick line in Figure 4(a) is the ensemble average of the 



field quadratures. This ensemble average, corresponding 
to an unobserved, dissipative evolution, remains close to 
zero, since, on average, the observer will detect the two 
symmetry-broken states an equal number of times. 

The second panel (b) in Figure 4 shows the von Neu- 
mann entropy S(p) = — Tr(/olog 2 p) for the oscillator 
state. In general, the entropy for the individual trajec- 
tories are decreasing towards zero, again consistent with 
the cavity approaching a pure state. The average (dis- 
sipative) evolution, however, does not have decreasing 
entropy. This is another signature of the fact that, for 
specific trajectories the observers knowledge is sufficient 
to distinguish the two symmetry-broken states. 

If the initial state was a pure superposition of the two 
many-body symmetry-broken states, the oscillator would 
be in an equal binary mixture of the two coherent states 
| a) and |— a) and we would expect 5=1. However, in 
the case shown here, 5(0) w 1.6 indicating that the initial 
state is not well-described by two orthogonal symmetry- 
broken states. This is not particularly surprising in this 
case, since the system is of moderate size (L = 32) and 
the chosen initial state is close to the first order phase 
transition shown in Figure 3. 

Note, that 5 is also a measure of the number of states 
in a Schmidt-decomposition of the state into oscillator- 
and spin-chain sub systems. Since the oscillator entropy 
is small for an accurate representation of the state (the 
ground state energy is accurate to 5 decimal places) we 
conclude that, in this case, the non-local interactions do 
not present a major problem for the MPS techniques. 

The third panel (c) in Figure 4 shows the expecta- 
tion value of the symmetry-operator P = n (g> a\ . . . ® . 
The figure clearly shows that the initial ground state 
is anti-symmetric and that the symmetry is lost in the 
subsequent time-evolution where one of the symmetry- 
broken states is selected. 



V. CONCLUSION 

In this paper we have applied matrix product state- 
methods to investigate ground-state and dynamic prop- 
erties of the Dicke-Ising model. The Dicke-Ising model 
not only has interactions between spins, but also a global 
cavity mode, interacting uniformly with all the spins. 
This implies that entanglement can exist over arbitrary 
distances in the system which, in turn, can be a chal- 
lenge for the numerical treatment of the problem. Wc 
have shown that the matrix product state method with a 
Krylov subspace expansion of the propagator permits ef- 
fective numerical treatment of the problem, and we have 
identified the phase transitions in the model predicted by 
mean-field theory. 

Near the second order phase-transition the mean-field 
solution for the average photon number is readily repro- 
duced, while the photon number variance deviates signif- 
icantly from the mean-field results. The increased vari- 
ance is ascribed to the formation of a superposition of 
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two coherent states near the phase transition in contrast 
to the even statistical mixture of coherent states pre- 
dicted by mean-field theory. Near the first order phase 
transition there are significant and non-trivial finite-size 
effects causing a weak thermal-like excitation of the cav- 
ity mode below the phase transition, which approaches 
the vacuum state only for very long spin chains. 

The interplay of quantum mechanical measurements 
and many-body physics is an interesting combination of 
thermodynamics and information theory and to investi- 
gate the phase transition dynamics further, we have sim- 
ulated the conditioned stochastic dynamics of the Dicke- 
ring system when it is subject to homodyne detection of 



the field leaking from the cavity. 

Quantum mechanical measurements serve as an at- 
tractive model for understanding the symmetry-breaking 
mechanism. Using this approach we can not only explain 
why the symmetry-breaking happens, but also the asso- 
ciated time-scale. As clearly shown by our simulations, 
the measurement back action causes a dynamic symme- 
try breaking, as the observer gradually distinguishes be- 
tween the cavity output from the two candidate (mean- 
field) symmetry-broken states. 

The authors acknowledge support from the EU inte- 
grated project AQUTE. 
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